Magnetic screening properties of superconductor-ferromagnet bilayers 
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We study theoretically the magnetic screening properties of thin, diffusive superconduc- 
tor/ferromagnet bilayers subject to a perpendicular magnetic field. We find that the effective pene- 
tration depth characterizing the magnetic response oscillates with the thickness of the ferromagnetic 
layer on the scale of the ferromagnetic coherence length. 
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While superconductor-normal metal (SN) struc- 
tures have been intensively studied for decades, 
superconductor-ferromagnet (SF) structures have only 
become accessible recently because of the much reduced 
length scales in ferromagnets. Due to their incompatible 
spin properties, the proximity effect between a singlet su- 
perconductor and a ferromagnet leads to a variety of un- 
usual phenomena 1,2 . Through the exchange field acting 
on the electron spins in the ferromagnet, Cooper pairs ac- 
quire a finite momentum 8k which leads to an oscillatory 
behavior of the anomalous Green function 3 -. Observable 
consequences are, e.g., a non-monotonic dependence of 
the transition temperature 4 *^ and the density of states 
at the Fermi level^ on the thickness of the F layer in 
SF bilayers, and the possibility of 7r-Josephson junctions 
at certain thicknesses of the F layer in superconductor- 
ferromagnet-superconductor (SFS) trilayer o 9 ' 10 ^ 11 . 

While most experiments on hybrid systems use resis- 
tive measurements, screening of an external magnetic 
field offers an alternative tool to study the proximity 
effect. These measurements probe deeply into the su- 
perconducting state because they provide both the mag- 
nitude and temperature dependence of the effective su- 
perfluid density. Various configurations for the magnetic 
response can be considered. The magnetization of SN 
hybrids with a magnetic field applied parallel to their in- 
terface has been addressed theoretically in Ref. 12 , with 
still debated experimental results in the case of SN cylin- 
der o 13 ' 14 . Alternatively, the screening properties of thin 
films can be probed by measuring the mutual inductance 
of two coils positioned on opposite sides of the sam- 
pl o 15 ' 16 . The mutual inductance can be related to the 
complex conductivity of the film which in turn can be 
related to the screening length A or the superfluid den- 
sity ps- To be precise, in SF bilayers, one measures the 
superfluid density ps <x A~ 2 integrated over the width of 
the bilayer or an effective screening length 

Kl = dt I * dx \-\x), (f) 

J-ds 

where ds and (If are the thicknesses of the supercon- 
ducting and ferromagnetic layer, respectively, and x is 
the coordinate normal to the interface. First experimen- 



tal results on Nb/Ni bilayers have been reported using 
this setup in Ref.—, where a non-monotonic dependence 
of the effective screening length A~^ on the thickness of 
the Ni layer has been observed. 

In this paper, we study the screening length A c ff of 
a SF bilayer subject to a weak perpendicular magnetic 
field. The main assumptions are that (i) the exchange 
field h in the ferromagnet is much larger than the super- 
conducting order parameter A, (it) the system is in the 
dirty limit and, thus, the Usadel equation 1 ^ can be used, 
(Hi) the screening length A er r is much larger than the 
thickness d = ds + (If of the bilayer, and (iv) the width 
ds of the superconducting (S) layer is smaller than the 
superconducting coherence length £s = \J D 5 / '(2ttT c0 ) , 
where Ds is the diffusion constant and T c q is the transi- 
tion temperature of the bare S layer. Our main result is 
that the screening length displays an oscillatory behavior 
with the thickness of the ferromagnet. 

Due to the normalization condition g = 1 of the quasi- 
classical Usadel Green function g, it can be parametrized 
by an angle 6 such that the normal Green function G — 
cos 6 whereas the anomalous Green function F = sin 9. 
The system is then described by four coupled equations 
in terms of the angles 63 on the S side of the SF interface, 
6*o on the ferromagnetic (F) side of the SF interface, and 
9p at the ferromagnet- vacuum interface. 

The Usadel equation of the F layer, —DpV 2 6 + 
2ihsm8 — 0, can be integrated to yield 

2Viy = / dff - (2) 

Je F y/COStlF — cost) 

where y — gJf/Cf and £p = \jDpjh is the ferromag- 
netic coherence length, with the diffusion constant Dp of 
the F layer. The boundary condition imposing current 
conservation at the SF interface 1 ^ can be expressed as 

sin(6» s - 6 ) = 2\fi(3y/co&0 F - cos6» , (3) 

where /3 = Rtap/^F. Here Rb is the interface resistance 
per square, and <jp is the conductivity of the F layer. In 
the limit ds -C £,s, the Usadel equation of the S layer, 
—Ds^7 2 + 2oj sin 8 = 2ACOS6 1 , where u> is a fermionic 
Matsubara frequency, can be simplified by an expansion 



2 



in small spatial variations of the angle 9 across the S layer 
combined with the boundary condition ©. One obtains 

u> sin 6 s + lyfia-yj cos Op — cos 9q — Acos#s, (4) 

where a = Ds&f / '(^sds^F) and as is the conductivity 
of the S layer. Finally, the self-consistency equation for 
the order parameter A reads 



A = ttTAbcs^ 



E 



sm tls 



(5) 



where Abcs is the BCS coupling constant. 

In diffusive superconductors, the screening length 
A describes the local (London) current response 2 ^ to 
a vector-potential, j = — 1/{^\ 2 )A, where A~ 2 = 
[2i:Tjj, as /h) F 2 is proportional to the superfluid 
density. Here /xq is the vacuum permeability. In SF bilay- 
ers, the effective screening length is related to the angles 
9 through the equation 
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/ dx sin 2 9(xj) 
o 



,(6) 



where 7 = <Tf£f /(csds)- Using the Usadel equation 
of the F layer, the integral over x can be traded for an 
integral over 8, namely dx = j^( cos 0f ~ cos 9)~ 1 / 2 d9, 
ranging from dp to 0o- 

In general the set of equations ©-© can be solved 
numerically only, but in some limiting cases an analytical 
solution is possible. At T = 0, a simplification occurs 
because the sums over w can be replaced by integrals, 
and subsequently the integration over w can be traded 
for an integration over #5 using the Usadel equation 21 . 
It is then sufficient to solve the Usadel equation at u> = 
for 6s(0). In particular, using Eqs. {2J and the 
Usadel equation of the S layer 0$ can be brought into 
the form (cu + F(9s))sm9s = A cos 83, yielding duj = 
-(A/ sin 2 8 S + F'(8 s ))d8 s . Using this trick, the zero 
temperature gap A is given as 
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? ) sin 9 s 



(7) 



where Ao is the zero-temperature gap of the bare S layer, 
a result which can then be used to compute A^ . 

In the following, we provide analytical results for the 
effective screening length in two limits, namely (i) for 
a system without barrier (3 = and (ii) for a system 
with a strong barrier (3 ^> 1. In both cases, solutions are 
presented for small parameters a. For convenience, we 
introduce the notation x = (1 + i)x for x = a, ft, y. 

In the absence of a barrier /3 = 0, the boundary condi- 
tion imposes that the angles on both sides of the SF 
interface, 8q and #5, are the same. 



If dp S> £f, the angle Of is small, and Eq. |(5J) yields 
dp = 8 tan Of- exp[— y\. Thus, we can simplify Eq. ^ to 
yield 



9 I tan 2 — 

10 sin 9 S + 2a sin 1 - 8e~ 2 ^ , „ , 4 



A cos S . (8) 



Treating a <C Ao perturbatively, one finds 6*5 (0) = § + 
80s, where S9 S = — v/2(a/A ) (l - 16e- 25 (3 - 2^2)), 



and 



SA = 



-2a(\/2-ln(l + \/2) 
16 



(9) 



%t[ae~ 2 y] (31n(l + y2) + 4-5%/2 



This solution describes gapless superconductivity with 
a finite density of states v(0) at the Fermi level in the 
superconductor that oscillates with the thickness of the 
ferromagnet: v(0) = — vq^R[S9s], where vq is the density 
of states in normal state. The equation for the screening 
length takes the form 

i = l-(a 1 -a 2 (cos2y + sin2 2 /)e- 2 f)^- (10) 

2 ^2 / _ 2 „ „ a 

+ — — 7 l + a 3 ye y cos2y-a 4 — 
3tt V A 

where Aq~ 2 = nfj, asA /h is the inverse screening length 
of the bare S layer at zero temperature, and the coeffi- 
cients dj are positive^ 2 -. Both, the contributions to the 
effective screening length from the S layer and from the 
F layer (ex 7), oscillate on the length scale of the ferro- 
magnetic coherence length. 

If on the other hand dp <C £f, the variation of the 
angle 8 is small across the F layer and, thus, 9s — 9 p <C 1. 
Using Eq. ([2|), one obtains 83 — 8f coshy. Inserting this 
relation into the Usadel equation of the S layer results in 



2iay + —ay cos 9s ) sin 6s — Acos0s. (11) 



We find S9 S = ~2iay/A Q and SA = -fay 3 . Note that 
because 9s (0) = §-H</>, where real, the density of states 
possesses a gap in this regime. Using Eq. (fTTj) to convert 
the integral over oj to an integral over 9s, the screening 
length is given by 
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7r ,„ 16 . a o 
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(12) 



Eq. (fTSj) predicts an increase in A~g as long as dp < 
before it starts to decrease. The regime dp ~ 
£f connecting the results Eq. (|10p and (fl"2"|) is treated 
numerically (see below). 

In the opposite limit of a strong barrier, /? > 1, both 
0o and 9f are small, if the F layer is not too thin, y ^> 
Eq. @ then yields 9 = 9f coshy, and, using the 



3 



boundary condition, the Usadel equation of the S layer density of states, v(0) — vqol/{(}A.q), while 
can be rewritten as 



a 

w + - 



- cos 9s I sin 9s = Acos#s. (13) 



P \/2i/? 2 tanhy 

Eq. (|13p yields 59s = —o-l (/?Ao), and thus no gap in the The screening length is given as 
I 



5A = ~ + ^[(l-i)cothy\. (14) 
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Again the effective screening length oscillates on the scale 
of £p. However, these oscillations are suppressed due to 
the large barrier. For y ^> 1, the oscillatory function in 
the contribution of the S layer has the same form as the 
one in (p~0|) whereas the oscillating part of the contribu- 
tion from the F layer is proportional to ye~ 2y sin(2j/). 

In the case of a very thin F layer, y <§C /3 _1 , the varia- 
tion of the angle 9 is small across the F layer, see above. 
The boundary condition at the SF interface then simpli- 
fies, and the Usadel equation of the S layer yields 

(u; + 2iay + 4a f3y 2 cos 9s) sin 9s — Acos9s- (16) 

We find 59s — — 2iay/A and 5A = —naPy 2 . As for 
the case without a barrier, the density of states in the 
thin film regime is gapped. Using Eq. (TT6")) to convert 
the integral over lo to an integral over 9s, the screening 
length is given by 
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The inverse screening length increases in the very narrow 
regime dp < £,p/(P^ 2 s )- 

Thus, we find oscillations of the screening length both 
in the absence of a barrier and in the presence of a strong 
barrier. The amplitude of oscillations in the latter case is 
suppressed, however. In both cases, analytic results can 
be found for small thicknesses y < y* and large thick- 
nesses y > y* , where y* ~ miri{l,/3 } denotes the po- 
sition of the first strong minimum. The vicinity of this 
minimum is not accessible to analytic solution. 

To find a solution in this regime, we note that Eq. ([2]) 
yields a general relation between 9q and 9p, namely 
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where cn is the Jacobi elliptic function. Using Eq. (TTSJ), 
the boundary condition ([3]) yields 9s as a function of 9f- 
Inserting this solution into the other equations, the set 
of coupled Eqs. and (O can then be solved numeri- 
cally. The thickness dependence of the low-temperature 




FIG. 1: Oscillation of the inverse screening length 1/Ae ff as 
a function of cLf at temperature T = 0.1T c o. Here we use the 
parameters a — 1.2, /3 = 1, and 7 = 0.6. The inset magnifies 
the weak maximum at d,F ~ 2.5£f- 



screening length is displayed in Fig. [T] The minimum at 
d-F ~ £f is clearly visible whereas further oscillations at 
larger dp are very small. 

Furthermore, the numeric solution allows one to de- 
scribe the temperature dependence of the screening 
length. Fig. [5] shows the finite temperature curves for 
different values of dp/^p. The oscillations of the zero- 
temperature screening length mirror the oscillations of 
the critical temperature as well as the slope of the screen- 
ing length close to T c . 

In the vicinity of the critical temperature T c an an- 
alytic solution is possible for all parameter values. For 
the SF bilayer, the critical temperature is given by the 
solution of the equation^ 



T c0 V2 



5R 



1 



(19) 



where the (complex) relaxation time r s reads 

_ 1 \f2~i ata.nh(\/2i y) 

1 + V2i/3tanh(^y)' 

Eqs. (|19l20p are obtained by linearizing ©-J!]) in small 



(20) 
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9 close to the transition. For r" 1 small, T c — T c q — 
jdt [t~ x ] ■ For \t~ x \ = Aq/2, the transition temperature 
vanishes according to Eq. (fl9|) . Note, however, that for 
large r s _1 the transition typically becomes first order—. 



Expansion of equations (0)-(O up to cubic order 
then yields the temperature dependence of the screen- 
ing length A e ff close to T c . Namely, 



hirT r 



A 2 (T)5i 
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2y + sinh(2y) 



4v2i (cosh y + (3 sinh y) 2 
I 



1 



2ttT c t s 



(21) 



We see that the contribution of the F layer to A~^ dis- 
plays an oscillatory dependence on its thickness. Fur- 
thermore, both T c and A(T) oscillate. In particular, 



A 2 (T) = 4^T c (l-?R[(27rr e T s )- 1 ^( 1 )(0)]) 
T c -T~ -S}[^)(.)-/(rJ,y)*(3)( 2 )] ' 

at T < T c , where z = \ + (2irT c T a )- 1 and 



(22) 



(23) 



l9 2y + sinh(2y) \ 

4(1 + 2/3tanhy) 2 - \Jl 
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Note that the simple relation A -2 oc A 2 does not hold in 
the presence of the F layer. The slope of A~^ close to T c 
has its own dependence on the thickness of the F layer 
and the relaxation rate r. . 
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FIG. 2: Temperature dependence of \~2 . The parameters 
used are the same as in Fig. and four different thicknesses 
are shown: <1f/£f = 0.5, 0.9, 1.3, and 3. The corresponding 
critical temperatures are T c (0.5) = 0.36T c0 , T c (0.9) = 0.17T c0 
(close to the minimum), T c (1.3) = 0.23T c0 , and T c (0.5) = 
0.31T c o (close to the asymptotic value for <1f/^f S> 1). 



In conclusion, we have shown that the screening length 
of SF bilayers displays a oscillatory dependence on the 
thickness of the F layer. Analytic solutions have been 
found in various regimes and a general solution has been 
determined numerically. The obtained non-monotonic 
dependence of the screening length has been observed 
experimentally^. Our method can be easily extended to 
more complicated situations such as multilayers where 
unusual features of the proximity effect have been pre- 
dictedii 2 -. 
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